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Abstract 

We study the blackbody spectrum from slabs of three-dimensional metallodielectric photonic 
crystals consisting of gold nanoparticles using an ab initio multiple-scattering method. The spectra 
are calculated for different photonic-crystal slab thicknesses, particle radii and hosting materials. 
We find in particular that such crystals exhibit a broadband emission spectrum above a specific 
cutoff frequency with emissivity of about 90%. The studied photonic crystals can be used as 
efficient selective emitters and can therefore find application in thermophotovoltaics and sensing. 
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The main feature of photonic crystals is the ability to tailor the photon density of states 
and this way control the spontaneous emission of light, aiming at the realization of new opto- 
electronic devices. In this context, there has been considerable effort to design and fabricate 
photonic crystals which allow for control of thermal emission of light, i.e. thermally driven 
spontaneous emission, promising applications in imaging, sensing and most importantly, in 
thermophotovoltaics (TPV) . 1 i 2 i 3 i 4 i 5 i 6 i 7 i 8 i 9 Control of thermal emission can also be achieved 
by means of microstructured engineering on silico n 10 ' 11 ' 12 ' 13 or metal surfaces.— Depending 
on the type of application, photonic crystals and structured surfaces can act as narrow- or 
wide-band, directional or isotropic thermal emitters. For example, in TPV applications^ a 
quasi-monochromatic emission is preferable whilst in radiation cooling 1 — a broad emission 
spectrum is desired. In this work we investigate the emission properties of three-dimensional 
metallodielectric photonic crystals consisting of gold nanoparticles. We find, in particular, 
that the emission spectrum of these crystals can be such that photons are emitted in all 
directions only when their energies lie above a specific cutoff frequency, with emissivity as 
large as 90%. 

Photonic crystals of spherical scatterers have been theoretically studied using multiple 
scattering theor y 17 ' 18 which is ideally suited for the calculation of the transmission, reflection 
and absorption coefficients of an electromagnetic (EM) wave incident on a composite slab 
consisting of a number of planes of non-overlapping particles with the same two-dimensional 
(2D) periodicity. For each plane of particles, the method calculates the full multipole expan- 
sion of the total multiply scattered wave field and deduces the corresponding transmission 
and reflection matrices in the plane-wave basis. The transmission and reflection matrices 
of the composite slab are evaluated from those of the constituent layers. Having calculated 
these matrices one can evaluate the transmittance T(u,9,(fi), reflectance 1Z(oj , 9 , <ft) , and 
from those two, the absorbance A(u>, 6, (f) of the composite slab as functions of the incident 
photon energy fijjj and incident angles 9 and <fi. Transmittance and reflectance are defined 
as the ratio of the transmitted, respectively the reflected, energy flux to the energy flux as- 
sociated with the incident wave. The method applies equally well to non-absorbing systems 
and to absorbing ones. In terms of speed, convergence and accuracy, the multiple scattering 
method is the best method to treat photonic structures of spherical particles. The emittance 
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S(u, 6, 0) is calculated indirectly by application of Kirchoff 's law , 19 i 20 i.e. from 

E(u, 6, 0) = A(u, 9, (j)) = 1 - K(w, 6, (p) - T(w, 6, 0) (1) 

Note that recent direct calculations of the thermal emission from photonic crystals have 
verified the validity of Kirchoff 's law for the case of photonic crystals.— 

In this work we deal with gold spheres with radius of a few nanometres. Traditionally, 
the gold spheres are treated as plasma spheres whose dielectric function is given by Drude's 
formula 

£ * M = 1 " CO 
where uo p stands for the bulk plasma frequency of the metal and r is the relaxation time 
of the conduction-band electrons. The dielectric function e(uj) = £\(oS) + i e 2 (uj) for bulk 
gold has been determined experimentally by Johnson and Christy^ 2 - over the range from 
hu = 0.64 eV to hw = 6.60 eV. Taking frr" 1 = Ht^ 1 = 0.027 eV and hw p = 8.99 eV-2^ 
Eq. (j2j) reproduces satisfactorily the experimentally determined 6i(o>) over the whole of the 
above frequency range, and the experimentally determined 62(00) at low frequencies. r& refers 
to the collision time of conduction-band electrons in the bulk metal; due to the finiteness of 
a nanosphere, the collision time r of an electron in a sphere of radius S is decreased due to 
scattering at the boundaries of the sphere. It is given by 



T- 1 



t^+vfS- 1 (3) 



where vf stands for the average velocity of an electron at the Fermi surface. For gold 
hvp = 0.903 eV-nm , 2 ?' 24 so that for a sphere of radius S — 5 nm, using frr^ 1 = 0.027 eV, 
we obtain fir~ x = 0.21 eV. Therefore, a more realistic description of the dielectric function 
of a gold sphere would b o 23 ' 24 



e s {u) = e{u) + - p — - p (4) 

uj{uj + 1 r b ) uj(uj + ir ') 

where e(u) is the experimentally determined dielectric function of bulk gold. We note that 
65(0?) as given above, is different from the approximation to the dielectric function of the 
gold sphere obtained from Eq. (j2j) with hio v = 8.99 eV and frr -1 = 0.21 eV only because the 
imaginary part of the dielectric function is not represented well by Eq. (J2J) for hw > 1.5 eV. 
Taking into account the finite size of a metal nanoparticle for the determination of the actual 
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dielectric function, i.e. Eq. (JJJ), has successfully reproduced experimentally obtained light 
absorption and scattering spectra of monolayers of gold nanoparticles.— 

To begin with, we consider an fee photonic crystal whose lattice sites are occupied by 
gold nanospheres of radius 5 = 5 nm. The lattice constant of the crystal is a = 19.11 nm 
corresponding to a volume filling fraction occupied by the spheres, / = 0.3. The spheres 
are supposed to be suspended in air. The crystal is viewed as a succession of planes of 
spheres (layers) parallel to the (001) surface of fee. In Fig. ^we show the transmittance, 
reflectance and absorbance vs energy for light incident normally on a finite slab of the crystal 
consisting of 128 layers of spheres. We observe that for energies above 2.2 eV almost all light 
(92%-98%) that is incident on the crystal slab is absorbed from the gold spheres. Therefore, 
from Kirchoff 's law, i.e. £(u, 6, <p) = A(u>, 6, <p), we can infer that the blackbody radiation is 
strong above this energy. Indeed, the blackbody radiation intensity of the photonic crystal 
is given by 



ks is the Boltzmann constant and c is the speed of light in vacuum. In a metal particle 
of nanoscale radius, there are two principal sources of light emission/ absorption: the first 
one is the dipole oscillations due to the surface plasmon resonances and the second is the 
interband transitions that take place in bulk gold contributing to its dielectric function. In 
terms of the first source of emission, the crystal can be viewed as a lattice of interacting 
oscillating dipoles where each one of the dipoles emits light around the surface plasma 
resonance frequency ujsp = w p / v^3- The interaction between the dipoles leads to the creation 
of a band of resonances around u>sp- In two-dimensional systems, this band of resonances is 
manifested as a strong absorption/ emission peak 28 whilst in three-dimensional systems is 
manifested as two absorption/ emission peaks (corresponding to two bands of resonances) 
with a reflectivity peak (stemming from a photonic band gap) lying between them.— In 
Fig. n we cannot clearly identify surface plasmon peaks as these must have been submerged 
into the absorption edge due to interband transitions. In any case, both contributions 
(plasmon+interband) lead to a strong absorption/ emission plateau above 2.2 eV. 

As the system considered here can potentially found application in TPV devices we are 




(5) 



where Ibb{u,9,4>,T) is the blackbody radiation intensity (Planck distribution) 




(6) 
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interested in the spectral hemispherical (SH) radiative properties of the photonic crystal. 
More specifically, we are interested in the SH emissivity s(u) which is the ratio of the SH 
emissive power of a photonic-crystal slab to the SH emissive power of a perfect blackbody 
at the same temperature T. For a slab infinitely extended in two dimensions: 



Note that £(u),9,<j)) in the above equation is the arithmetical average of both polarization 
modes. In Fig. |2 we show the SH emissivity from slabs of the photonic crystal of Fig. ^ 
for different slab thicknesses (number of layers). It is evident that as the number of layers 
increases the SH emissivity increases accordingly until it reaches a saturation plateau for 
slabs of 64 layers and above. Indeed, for energies above 2.2 eV the emissivity curves for 64 
and 128 layers are practically the same. This fact implies that, if we measure the emission 
from one of the surfaces of a 128 layers-thick slab, then the emitted radiation must be coming 
from the 64 outmost layers (relative to a given surface) whilst radiation coming from the 
64 innermost layers is almost totally absorbed from the 64 outmost layers by the time it 
reaches the surface of the slab. The most important finding of Fig. El is the fact that the 
photonic crystal behaves more or less as gray body (emissivity ranging from 88% to 92%) for 
energies above 2.2 eV while, at the same time, emits small amounts of radiation for energies 
below this cutoff. A similar emission spectrum is observed for the Salisbury screen^ but 
with significantly lower values of the SH emissivity. 

So far, we have assumed that the gold nanospheres are suspended in air. In order to 
provide a manufacturable structure it is necessary to examine the case where the spheres 
are embedded in a host material of dielectric function e^. In Fig.Elwe show the SH emissivity 
from a 64-layer slab of an fee photonic crystal of gold nanospheres with / = 0.3 where the 
spheres are surrounded by air (e^ = 1), silica (e^ = 1.88) and gelatine (e^ = 2.37). The 
introduction of a host material does not change the picture drastically except that it lowers 
the emission cutoff energy. This is due to the lowering of the surface plasmon frequency 
according to ojsp = 0J P / \/l + 2e^ when a nanosphere is placed in a host medium e/j, which 
results in strong thermal emission at lower energies. As a result of the energy shift of the 
surface plasmon resonance, the corresponding emission peak moves away from the interband 



J Q 27r d(j) f^ 2 dOIpcioj, 9, (ft, T) cos 9 sin 9 
J 2w d(j) J^ 2 d9I BB (ui, T) cos 9 sin 9 
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emission region and a distinct local minimum appears. This is evident from the emissivity 
spectra of Fig. |3] for silica and gelatine host materials. The shift of the surface plasmon 
energy with respect to the host material allows for tuning of the cutoff energy so that it 
coincides with the band gap of the photodiode of a TPV device. 

Next we study the effect of the volume filling fraction / on the emission properties of the 
photonic structure under study. In Fig. 0] we show the SH emissivity spectra for photonic 
crystals of different values of /. We have kept the lattice constant the same and changed 
the sphere radii accordingly. Note that for each value of / we obtain a different value of 
hr^ 1 from Eq. (j^J. We observe that the maximum emissivity in the region above the cutoff 
energy is achieved for / = 0.3. Since our calculations show that the SH reflectivity assumes 
higher values for / = 0.5,0.7, the radiation emitted from the innermost layers is reflected 
back and reabsorbed before it reaches the surface of the slab. 

Finally we address the issue of spatial order/ disorder of the photonic crystal under 
investigation. As it has been both theoreticall y 30131 and experimentally 32 shown, the pres- 
ence of disorder does not changes practically the absorption/ emission properties of three- 
dimensional photonic crystals of metal nanoparticles since these quantities mostly depend 
on single sphere properties such as surface plasmon resonances, and less on the particular 
spatial arrangement of the spheres. Note that this not true for one- or two-dimensional 
structures, e.g. a linear chain or a monolayer of nanospheres, where the presence of dis- 
order changes dramatically the absorption/ emission spectra.— So, since we deal with a 
three-dimensional structure, a disordered crystal might be easier and cheaper to fabricate. 

In summary, we have shown that a metallodielectric crystal consisting of gold nanospheres 
can act as a 90% gray body for energies above a cutoff. The cutoff frequency can be tuned by 
appropriate choice of the material hosting the nanospheres. Optimal emissivity is obtained 
for moderate volume filling fractions and slab thicknesses. 

ACKNOWLEDGEMENTS 

This work was supported by the 'Karatheodory' research fund of University of Patras. 



Electronic address: vyannop@upatras.gr 



1 C. M. Cornelius and J. P. Dowling, Phys. Rev. A 59, 4736 (1999). 

2 S. Y. Lin, J. G. Fleming, E. Chow, J. Bur, K. K. Choi, and A. Goldberg, Phys. Rev. B 62, 
R2243 (2000). 

3 J. G. Fleming, S. Y. Lin, I. El-Kady, R. Biswas, and K. M. Ho, Nature (London) 417, 52 (2002). 

4 M. U. Pralle, N. Moelders, M.P. McNeal, I. Puscasu, A. C. Greenwald, J. T. Daly, E. A. Johnson, 
T. George, D. S. Choi, I. El-Kady, and R. Biswas, Appl. Phys. Lett. 81, 4685 (2002). 

5 S. Y. Lin, J. G. Fleming, and I. El-Kady, Opt. Lett. 28, 1909 (2003). 

6 I. Celanovic, F. O' Sullivan, M. Ilak, J. Kassakian, and D. Perreault, Opt. Lett. 29, 863 (2004). 

7 A. Narayanaswamy, and G. Chen, Phys. Rev. B 70, 125101 (2004). 

8 S. Enoch, J. -J. Simon, L. Escoubas, Z. Elalmy, F. Lemarquis, P. Torchio, and G. Albrand, 
Appl. Phys. Lett. 86, 261101 (2005). 

9 M. Florescu, H. Lee, A. J. Stimpson, and J. Dowling, Phys. Rev. A 72, 033821 (2005). 

10 H. Sai, H. Yugami, Y. Akiyama, Y. Kanamori, and K. Hane, J. Opt. Soc. Am. A 18, 1471 
(2001). 

11 S. Maruyama, T. Kashiwa, H. Yugami, and M. Esashi, Appl. Phys. Lett. 79, 1393 (2001). 

12 F. Marquier, K. Joulain, J. P. Mulet, R. Carminati, and J. J. Greffet, Opt. Commun. 237, 379 
(2004). 

13 M. Laroche, F. Marquier, R. Carminati, and J. J. Greffet, Opt. Commun. 250, 316 (2005). 

14 F. Kusunoki, T. Kohama, T. Hiroshima, S. Fukumoto, J. Takahara, and T. Kobayashi, 
Jpn. J. Appl. Phys. 43, 5253 (2004). 

15 M. Zenker, A. Heinzel, G. Stollwerck, J. Ferber, and J. Luther, IEEE Trans. Electron. Devices 
48, 367 (2001). 

16 C. W. Hoyt, M. P. Hasselbeck, M. Sheik-Bahae, R. I. Epstein, S. Greenfield, J. Thiede, J. Distel, 
and J. Valencia, J. Opt. Soc. Am. B 20, 1066 (2003). 

17 N. Stefanou, V. Karathanos and A. Modinos, J. Phys.: Condens. Matter 4, 7389 (1992). 

18 N. Stefanou, V. Yannopapas, and A. Modinos, Comput. Phys. Commun. 113, 49 (1998); ibid 
132, 189 (2000). 

19 S. M. Rytov, Yu. A. Kravtsov, V. I. Tatarskii, Principles of statistical radiophysics, Vol. 3: 
Elements of random fields (Springer, Berlin, 1989). 

20 J. J. Greffet and M. Nieto-Vesperinas, J. Opt. Soc. Am. A 15, 2735 (1998). 

21 C. Luo, A. Narayanaswamy, G. Chen and J. D. Joannopoulos, Phys. Rev. Lett. 93, 213905 

7 



(2004). 

R. B. Johnson and R. W. Christy, Phys. Rev. B 6, 4370 (1972). 

S. Norrman, T. Andersson, C. G. Granqvist, and O. Hunderi, Phys. Rev. B 18, 674 (1978). 
F. Abeles, Y. Borensztein, and T. Lopez-Rios, Festkorperprobleme (Advances in Solid State 
Physics) vol 24 (Braunschweig: Vieweg) p. 93(1984). 

M. M. Wind, P. A. Bobbert, J. Vlieger, and D. Bedeaux, Physica A 143, 164 (1987). 

P. A. Bobbert and J. Vlieger, Physica A 147, 115 (1987). 

N. Stefanou and A. Modinos, J. Phys.: Condens. Matter 3, 8135 (1991). 

N. Stefanou and A. Modinos, J. Phys.: Condens. Matter 3, 8149 (1991). 

V. Yannopapas, A. Modinos, and N. Stefanou, Phys. Rev. B 60, 5359 (1999). 

A. Modinos, V. Yannopapas, and N. Stefanou, Phys. Rev. B 61, 8099 (2000). 

V. Yannopapas, A. Modinos, and N. Stefanou, Opt. Quant. Elec. 34, 227 (2002). 

K. P. Velikov, W. L. Vos, A. Moroz, and A. van Blaaderen, Phys. Rev. B 69, 075108 (2004). 



1.0 




Photon Energy (eV) 

FIG. 1: Transmittance, reflectance, and absorbance of light incident normally on a 128-layers thick 
slab of a fee crystal consisting of gold spheres (S = 5 nm) in air (e^ = 1), with / = 0.3. 
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FIG. 2: SH emissivity for different numbers of layers (slab thicknesses) of the photonic crystal 
described in Fig. ^ 
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FIG. 3: SH emissivity of a 64-layer thick slab of an fee crystal of gold spheres (S = 5 nm) in air 
(eft = 1 - solid line), silica (e^ = 1.88 - dashed line) and gelatine (e^ = 2.37 - dotted line), with 
/ = 0.3. We have assumed that the host medium surrounding the spheres also covers the whole 
space. 
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FIG. 4: SH emissivity of a 64-layer thick slab of an fee crystal of gold spheres (S = 5 nm) in air, 
for different volume filling fractions / (fixed lattice constant and different sphere radii). 
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